Influence of internal structure on the motion of test bodies in extreme mass ratio 

situations 
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We investigate the motion of test bodies with internal structure in General Relativity. With 
the help of a multipolar approximation method for extended test bodies we derive the equations 
of motion up to the quadrupolar order. The motion of pole-dipole and quadrupole test bodies is 
studied in the context of the Kerr geometry. For an explicit quadrupole model, which includes spin 
and tidal interactions, the motion in the equatorial plane is characterized by an effective potential 
and by the binding energy. We compare our findings to recent results for the conservative part of 
the self-force of bodies in extreme mass ratio situations. Possible implications for gravitational wave 
physics are outlined. 
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INTRODUCTION 



Approximate analytic treatments of the binary prob- 
lem in General Relativity are possible in certain lim- 
ited regimes. These include the weak field (post- 
Minkowskian) , weak- field and slow motion (post- 
Newtonian), and extreme mass ratio approximation. An 
important astrophysical example of the latter approxima- 
tion is a star inspiraling into a very massive black hole, 
i.e., an extreme mass ratio inspiral. The zeroth order ap- 
proximation in the extreme mass ratio case describes the 
motion of the lighter mass by a geodesic in the fixed back- 
ground spacetime generated by the heavier mass. An im- 
portant first correction derives from a perturbation of the 
background due to the small mass object, leading to so- 
called self-force corrections to the geodesic motion within 
the background metric. Recently this self-force program 
made progress, for example evolutions over thousands of 
orbits for a Schwarzschild background succeeded [l|, see 
also 0, Q . Also self- force corrections to the innermost 
stable circular orbit were derived whose frequency is 
an important observable of gravitational wave astronomy. 

Besides self-force corrections to the geodesic motion 
within a background metric, corrections due to the small 
object's spin and higher multipoles - encoding its finite 
size - arise. Multipolar corrections to the motion of test 
bodies were first derived by Mathisson [H, Q and Papa- 
petrou 0, and since then have been studied in the con- 
text of different approximation schemes by a number of 
authors (cj-il8j. It is important to systematically study 
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these corrections in astrophysically realistic situations, in 
particular as the self-force calculations are approaching 
the second order [2(| . In the present paper we study a 
truncation of the equations of motion at the quadrupole 
level, to which an independent derivation was given in our 
previous publication [21| with the help of the multipolar 
approximation scheme by Tulczyjew [8(. Our analysis in 
|2l| made clear that there is a considerable freedom to 
close the equations of motion at the quadrupole level. 
Besides imposing a supplementary condition on the spin, 
corresponding to a choice for the center of the object, one 
must devise a model for the quadrupole dynamics. 

In contrast to previous works [22-24], in the present 
paper we focus on astrophysically realistic models for 
the quadrupole, which derive from recent work on effec- 
tive actions for extended objects. This not only includes 
quadratic-in-spin corrections [25l [26j . but also tidal in- 
teractions [13, HH, see also [29j . In fact, the interaction 
of a black hole's tidal field with the quadrupole of an in- 
spiralling star may become very strong, as is indicated 
by the appearance of tidal disruptions. Quadrupole tidal 
interactions can be of more importance for gravitational 
wave astronomy in large mass ratio situations than spin 
effects - for the measurability of spin effects see f30j | - 

The present paper aims at extending the effective po- 
tential for test bodies with spin in a Kerr background [31] 
to include corrections from the mentioned quadrupole 
models (see also [H, |33| for alternative derivations, (34J 
for a char ged spinning test body in Kerr-Newman space- 
time, and [35|,[36| for Hamiltonian approaches). To carry 
out this program, it is necessary to identify various con- 
served quantities of quadrupolar test bodies. This is 
not straightforward and some of the emerging difficulties 
were already pointed out in our previous work [2l| . How- 
ever, with an effective action as a basis one can immedi- 
ately construct conserved mass and spin length parame- 
ters. Conserved quantities - associated with symmetries 
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(Killing vectors) of the background - were already found 
in [22( 1 and are actually the same as in the spinning test 
body (or pole-dipole) case. In a Kerr background this 
gives rise to a conserved binding energy and a total an- 
gular momentum of the test body. 

For aligned spins and circular orbits in the equatorial 
plane we compute the gauge-invariant relation between 
binding energy and total angular momentum. This com- 
plements recent results on this relation from full numeri- 
cal simulations for comparable masses [37j and from the 
conservative part of the self-force for extreme mass ra- 
tios [38] . While both of these results are valid for non- 
spinning binaries only, we include spin and quadrupole ef- 
fects (including tidal deformations) due to the smaller ob- 
ject in the extreme mass ratio case here. Furthermore, we 
compare our result with various post-Newtonian Hamil- 
tonians or potentials. 

The structure of the paper is as follows. In section 
HI! we briefly recapitulate some of our findings from [2l[ . 
We then rewrite the equations of motion for a gener- 
alized supplementary condition. Subsequently, in sec- 
tion IIII1 we work out the effective potential for spinning 
test bodies on equatorial orbits with aligned spin in a 
Kerr background. In section HVl we introduce our explicit 
quadrupole model. In particular we identify a mass-like 
quantity at this multipole order which is conserved in 
an approximative sense. Furthermore, we work out the 
effective potential in terms of a set of dimensionless pa- 
rameters, and discuss the values of these parameters for 
astrophysically relevant situations. This is followed by a 
discussion of the binding energy in section |Vj in which 
we also compare our results to corresponding ones in a 
self-force and post-Newtonian context. We draw our final 
conclusion in section |VT1 In appendix lAl we collect some 
misprints which we found in the literature regarding the 
effective potential of pole-dipole test bodies. Appendices 
IB] [El contain supplementary material regarding our cal- 
culations. In appendix [F] we summarize our notation 
and conventions and provide a brief overview of different 
quantities and units used throughout the work. 



II. BASIC EQUATIONS 
A. Equations of motion 

The equations of motion of an extended test body up to 
the quadrupolar order are given, see [HI for a derivation 
in the context of Tulczyjew's multipolar approximation 
method, by the following set of equations: 



5p a 
ds 
6S ab 



1 



-RabcdU S 



b Qcd 



V a RbcdeJ' 



bcde 



— = 2pK b ^ - -R [a cdE J b]cde 
ds 3 



(1) 
(2) 



Here u a := dY a /ds denotes the 4- velocity of the body 
along its world line (normalized to u a u a — 1), p a the mo- 
mentum, S — S ba the spin, and J abcd the quadrupole 



moment with the following symmetries: 



j[a 



bc]d 



o & j a - 



jabcd j^[ab][cd] jedab fo\ 

jcabd = Q ^ 



bed i jbcad 



Thus, J abcd has the same (algebraic) symmetries as the 
Riemann tensor. The corresponding stress-energy tensor 
of the test body can be written in the following (singular) 
form: 



ds 



- \Rede {a J h)edC 5 {i) 



-V C (S C K%)) - -Vd^e(J diab)C S (i) ) 



(5) 



From equation ^ it follows that the momentum is given 
by: 

<W ab 4 

p " = mv ? + Ub + u b R [a edeJ b]cde , (6) 

ds 3 
where we used m := p a u a . 



B. Conserved quantities 

We encountered already in [2l| , that for the pole-dipole 
case (J abcd = 0), the quantity 



El = PaC + ls ab V a C 



(7) 



is conserved if £ a is a Killing-vector, V( b £ a ) = 0. In 
[22L p. 2 10] it was further shown that this is a conserved 
quantity even at all higher multipole orders. 

Other conserved quantities depend on the supplemen- 
tary condition. For the spin length S given by 2S 2 := 
S a bS ab one obtains: 



s dS _ ldtP _ 1 SS ab 

ds 2 ds 2 ab ds 
2 

= S ab p a U b --S ab R a cdeJ bCde 



(8) 



It is easy to see that in the pole-dipole case the spin 
length S is conserved for the two well-known supplemen- 
tary conditions of Tulczyjew (p a S ab = 0), and Frenkel 
(u a S ab = 0). 

For the mass m, defined by m 2 := p a p a , the following 
relations hold: 



m = m 

dm _ 8p a p b 
ds 



X cab a 
2 + ~^PaUb + - Pa U b R^edeJ b]cde , (9) 

ds 3 



ds mm 



^_ + t R [a Jb]ede 

ds 3 



m SR b cde_jbcde 



6m ds 



(10) 



The last relation follows from an insertion of @ into the 



expression ^j^Pb^4 — and use of 



6p n 



am TT 

m-j=. Hence 
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in the pole-dipole case the mass to is conserved if one 
chooses Tulczyjew's spin supplementary condition. How- 
ever, for the Frenkel condition the mass m is conserved 
in the pole-dipole case. 

An extension of conserved spin length and masses to 
the quadrupole case is one of the main obstacles in the 
present work. We will come back to this in sec. II VI in the 
context of an explicit quadrupole model. 



C. Spin supplementary condition and equations of 
motion 



A spin supplementary condition (SSC) has to be im- 
posed to close the system of equations CO)-© even in 
the pole-dipole approximation. Here we rewrite the equa- 
tions of motion with the help of the following supplemen- 
tary condition: 



where 



S ab f b = 0. 



(*) 



(11) 



This condition allows for an easy transition between the 
widely used Frenkel (/„ = u a ) and Tulczyjew (f a = p a ) 
conditions. 

Taking the total derivative of (ITU) leads to a relation 
between p a and u a : 



1 



p f fi I 



*,bf „fl i gab h ^ f r>[a jb]cde 
U JbP + o -: nib**- cdeJ 

as 3 



(12) 



We now have the option to eliminate u a or p a from the 
equations of motion. 

Eliminating u a Insertion of (|12|) into the equation of 
motion for the spin ^ yields: 

ds ~ pff f ds 3 W f WgK cdeJ ■ [ 6) 

Here we defined the projector W b a with respect to the 
vector entering the supplementary condition i.e. 



w b a ■■= - K^, f a w b a = o, 

V s If 

W b a p b = 0, w?s cb = s ab . 



(14) 
(15) 



Now the spin dynamics is fixed, provided that is 
meaningful and J abcd is somehow given. This spin equa- 
tion of motion also manifestly preserves the spin supple- 
mentary condition. 

When u a is replaced in favor of p a , the following re- 
lation between the velocity and the momentum is often 
more useful than (TTZI): 



m- 



1 5 

m 2 ds 



S ab Pb ) 



3m 2 



R [a cdeJ h]cde 



Pb 



1 



6m 2 
to „ 



S ab V b RcdefJ cdef , 



III ~ 



-v 



S^Spb 
to 2 ds 



6m 2 



S ab V b RcdefJ Cdef 



3 J*/, 



4 ^ cde J b ^f b ' ^ SbC5fc 



■p f j 



pfff ds 



(17) 



(18) 



In the last equation we introduced p a := S b — p a p b /m 2 . 
The derivation of relation (|16l) at the quadrupole order 
is analogous to the one at the pole-dipole order given in 
(39} . Note that (fTB]) is only valid, if the supplementary 
condition (fTTj) holds. 

The mass quantity m := p a u a must be obtained from 
(|16p and u a u a = 1, as s is the proper time in our case. 
But one may as well use a different convention for the 
parameter s, namely u a U a = 1 where U a :— p a /m. Then 
it holds to = to, which makes the relation between u a and 
p a fully explicit, see also (22|. Now one can insert ([To]) 
into (UJ , which finally eliminates u a from the equations of 
motion. If desired, one can further decompose ([JJ into an 
equation for U a and to, but we will not explicitly follow 
this approach here. 

Eliminating p a Equation (|12p can be written as: 



1 



p b f b U a - S abS -^ + ^fbR^cdeJ^ 

ds 3 



(19) 



Insertion into the equation of motion for the spin ([2]) 
yields: 



~X [ ?X$Rf cde J<> cde , 



ds ff ds 3 ' 

where the projector X b is now given by 
u a fb 



va ra 

^b ■— °b ~~ f, 1 



faXg = 0, 



Xgu b = 0, x*s cb = s ab 



(20) 

(21) 
(22) 



As in the previous case this spin equation of motion man- 
ifestly preserves the spin supplementary condition. 
With the help of the relations 1 

faP a = mf a u a + S ab u a s -^ + y a u b R\ a cd eJ b]cde , (23) 

p a = mU «-^^ + §p»j!L-RV> cde J9]cde, (24) 

we can orthogonally decompose (UJ with respect to u a 



u a =p a + 



O oac ode r> 
Zb D tide 



Am 2 + S^S-fR cde / 



1 Here p% := 5£ — u a Uf, denotes the projector with respect to the 
velocity as usual. 
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into an equation for u a , and into an equation for m: 



Su a _ 1 
1 ds ~ 2 



RabcdU b S cd + p a 
4_ g fh 



6_ 

9 d~S 



v P 9 b S bc 5f c 



i f ff ds 



R [b cdeJ h]cde 



dm 
ds 



3 Pb uff f 

5u b S bc Sf c , Su b 4: f, 



-plVfR bcde J bode 



ds uf ff ds ds 3 uf ff 



9 R [b cde J 9]cde 



1 SRbcde jbcde 

6 ds 



(25) 



(26) 



Obviously u a u a = 1 is preserved. Observe that we have 
eliminated p a from all equations of motion. 



III. SPINNING TEST BODIES IN A KERR 
BACKGROUND 

In the following, we are going to study test bodies en- 
dowed with spin in the field of a rotating source described 
by the Kerr metric. 



A. Kerr metric 

In Boyer-Lindquist coordinates (t,r,9,<fr), the Kerr 
metric takes the form: 



ds 2 = 1 - 



2Mr\ , 2 

— r 



p z A 



, , . , , , 2 2a 2 Mrsm 2 6\ , , 
-p dff — sm ff I r + a + 1 



P 



(27) 



where M is the mass, a the Kerr parameter, and 

A = r 2 - 2Afr + a 2 , p 2 = r 2 + a 2 cos 2 6>. (28) 

The Kerr metric allows for two Killing vector fields given 
by: 

e = S?, T = 5%. (29) 

E J 

Furthermore, we have 

y/^g := V-det (g a b) = p 2 sin6. (30) 

B. Equatorial orbits for aligned spin 

In the following we are going to focus on equatorial 
orbits, i.e. 



and the case of aligned spin, with respect to the rotating 
background source, defined by 



S a0 = 0. 



(32) 



These assumptions leave us with the six remaining com- 
ponents {p t ,p r ,p' )> ,S tr ,S r 't > ,S'l' t }, which can be deter- 
mined from the following set of equations 



S ab Pb 



E = PaC 

E 

-J = PaC 



0; $abS 

l 

2 
1 



2S 2 , 

S a bW a C, 



PaP" 



(33) 



in terms of the mass m, the spin length S, the energy 
E, and the angular momentum J. Remember the cor- 
responding Komar quantities of the Kerr background at 
this point. The one belonging to the Killing vector £ a is 

E 

just the black hole mass (or energy), the other one, be- 
longing to £ a , yields the angular momentum. This moti- 

j 

vates to call E and J the energy and angular momentum 
of the test body, as they are based on the same isomc- 
tries, or Killing vectors. Note that from here on we use 
the Tulczyjew condition, i.e. the first equation in (|33[) . as 
supplementary condition in our multipole formalism. For 
this condition, we define the spin vector S a as follows: 



s a = ^-v abcd p b s cd = 

1 m im\/—g 



ibcd„ 



O = — 7] 

m 



abed 



PcSd 



i 



-<] 



PbScd, 

abed q 
Pc^d- 



(34) 



From this definition and the assumption in (|31j) it be- 
comes clear, that S a has only one non- vanishing compo- 
nent, i.e. 



S a = S ti S l 



e X a 



(35) 



Using the relation 2S a S a = -S a bS ab together with the 
definition of the spin length from (1331) , we obtain 

(36) 



- S e = ±S/^-g es . 



Notice that for usual spherical coordinates do points 
in the opposite direction as d z in the equatorial plane. 
Therefore S e < corresponds to a spin aligned to d z . 
Further discussion of the spin orientation is given in ap- 
pendix |51 As usual we absorb the sign by allowing for 
negative spin length S in the following. This in turn al- 
lows us to express the components of the spin tensor in 
terms of the spin length, i.e., 



S rt = - 



S^/—gee 



p<p 



SPcj, 



Q .pt _ S^J-gee _ Sp r 

— , Pr — i 

m^J—g m r 
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where in the last step we took into account that we are 
in the equatorial plane. 

The definitions of E and J in (|33|) allow us to express 
Pt and in terms of the constants of motion, i.e. 



E 



Pt 



P,b 



MS 
mr 3 



(J - aE) 



-J- 



aMS 



- J 



1 



MS 2 



(38) 
(39) 



With the last two expressions and the definition of the 
mass m from (|33[) we are able to express the radial com- 
ponent of the moment in terms of the constants of mo- 
tion, the test body mass, the test body spin, and the 
parameters characterizing the spacetime, i.e. 

(p r ) 2 = A (A 1 J 2 + A 2 JE + A 3 E 2 +A 4 ). (40) 

Here the functions A — A (M,m, S,r), Ai,...^ = 
{M,m, S, r, a) are explicitly given by: 

m 2 M 



\2 ' 



(S 2 M-m 2 r 3 y 
A x = S 2 M + 2mrSa + 2r 3 m 2 - 
2mr 4 S 



m 2 r 4 
M ' 



6r 3 Sm - Aam 2 r 3 - 2S 2 aM 



M 

-2S 2 ar — AmrSa 2 , 



(41) 
(42) 

(43) 



A 3 = 2a 2 m 2 r 3 



A 4 = 2r 3 S 2 m 2 



_ 2a 3 rSm + 2S 2 r 3 + 6r 3 Sam 

2S 2 ra 2 + S 2 Ma 2 

S 4 M - 
a 2 S 4 M 



M 



S 2 r 4 


2 4 

arm 


M 


+ M 


6 4 

r m 


2M 2 S 4 


M + 


r 



(44) 



AMr 2 S 2 m 2 



-2r 5 m 4 



2a 2 S 2 m 2 r 



? 4 4 

a m r 



M 



An alternative form of p r is the following one: 



E 2 



A 



a 2 




M 




r 


I- 


S x 


2 

)■ 




mr , 




a 2 




M 



no JE J 2 



(45) 



(46) 



r J mr 2 



7 
S 



A 








a 

- 4 


MS 


r 


mr 2 



. + r 



S 



m r 

\ 2 



A 



a= 1- 



A 



(47) 

(48) 
(49) 
(50) 
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(a) Parameter choice a = 0.5M, S = mM, J = 4m M. 
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(b) Parameter choice a = 0.5M, S = —mM, J = 4m Af . 

FIG. 1: Energy E = U±, in units of m, as a function of the 
radius r, in units of M. 



C. Effective potential (pole-dipole case) 

We may define the functions U + and C/_ by 



(ff=*(E-U + )(E-U-), (51) 
a 



so that 



U± =± i + M 2 $ 2 - ai J 2 + M^J_ 
m \/ a r 2 a 2 m 2 M 2 r a mM 



For p r to be a real number we need to have both E <U+ 
and E < U- , or both E >U+ and E > U- (under the 
assumption 2 that a > 0. For usual cases (and positive 
energy E > 0) the important relation is > [/+. This 
justifies to call U+ effective potential: The object can 
only move in the region where E > U+ and the turning 



2 We numerically checked that a > in the regime M < r < oo, 
-M < a < M, -Mm < 5 < Mm. 
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points are given by E = J7+, because then p r = - which 
implies u r — 0, see appendix [Ul Therefore the minimum 
ofU+, 



dU+ 
dr 



0, 



(53) 



together with E — U+ defines circular orbits. Thus, one 
can easily find circular orbit solutions by solving (|53|) for 
the radial coordinate r in terms of constants of motion 
and the other constant parameters. For this class of so- 
lutions the energy is constrained by E = U+. 

We have plotted the solutions for the energy E = U± 
for different spin configurations in fig. [TJ these can be 
directly compared to the ones found in [33} . 



IV. EXPLICIT QUADRUPOLE MODEL 

We are interested in the models for quadrupole defor- 
mations induced by spin given by [25|, eqs. (1) and (16)] 
and for adiabatic tidal quadrupole deformations given by 
[27l (19)], or by (2|| (5)]. These models are an extension 
of the point-mass action, see appendix iDl and the corre- 
sponding equations of motion have been worked out in 
a general fashion already in [40]. In particular, equation 
|40T (19)] provides a formula for the quadrupole J abcd . 
Other recent treatments (not based on an action princi- 
ple) of tidal effects in the context of General Relativity 
can be found in [4ll - [43j ] . Inspired by the abovementioned 
models we now consider 



jabcd 

where 



m 

77! 



l p [agfc]cd + L p [dQc]ba + 3 [a,Qb][c p d] 

m m mr 

(54) 



Q ab 

Qbcd 

E a b 
B a b 



= c ES 2S a e S be - fi 2 E ab , 

- _^V C v e B ba 

— '/ eaP E> ; 

m 

= -^RacbdP C p d , 



2m 2 



VaecdRbf Cd p e p f 



The quantities 0^52, /i 2 , and a 2 are assumed to be con- 
stants, and parameterize quadrupole deformations in- 
duced by the spin and by tidal forces of the spacetimc. 
Furthermore, E ab represents the gravito-electric tidal 
field, and B ab the gravito-magnetic (frame-dragging) 
tidal field, see, e.g., [44|. (The convention for B ab used 
in 28, 29] differs from the most common one adopted in 
the present paper by a factor of two.) Notice that the 
overall factor of m/m in J abcd makes the equations of 
motion rcparamctrization invariant. We choose the spin 
supplementary condition of Tulczyjew, 



S ab p b = 0, 



(55) 



as this condition is most convenient for the derivation of 
an effective potential (similar to the pole-dipole case). 



A. Spin length 



For the explicit model given by ([54]) we have 

2Q dc a Q dcb ) 



c c /no r\bde 

S-r- — ~ 7 — &ah\H deH 
US 0(7 2 



+fi 2 S ab E a d E bd + c ES 2E ad S ab S be S ed , (56) 

where (J5]) and the Tulczyjew condition were used. Fur- 
thermore, one has a symmetry of the combinations 



jpab jpa j^bd r\a r^bde / 1 ar\dcb 
& i & d& , W deW 7 Vdc V 1 

and the antisymmetry of 

S a b, S a dS S e b, 



(57) 



(58) 



under exchange of a and b. It immediately follows that 
the spin length is conserved, 



dS 
ds 



= 0. 



(59) 



This conservation law could be expected, as the action 
given in appendix |D]- which served as an inspiration for 
the present quadrupole model - possesses a symmetry 
under rotations of the body-fixed frame. But the action 
is only consistent to a certain power in spin, whereas the 
conservation found here made no reference to such an 
approximation. 



B. Mass-like quantity 

How a possibly conserved mass-like quantity is related 
to the usual masses m, or m, crucially depends on an 
explicit model for the quadrupole. In the present paper 
a conserved mass-like quantity is related to a parameter 
within the underlying action, which is a constant by as- 
sumption. However, the action is only consistent in an 
approximate sense, which translates here to the fact that 
we were only able to find a mass-like quantity which is 
approximately constant. 

We therefore introduce a multipole counting scheme 
by a symbol e of the form 



gab 



O(e ), 
0{e% 



O(e ), 



8p a 
ds 



SS ab 
ds 



0(e 2 ), 



= Otf), (60) 

jabcd = q ( £ 2)( 61 ) 



We will neglect O (e 3 ) terms in the following. Under these 
conditions, (fTU|) simplifies to 



d m 
ds 



1 SR a bcd_ jabcd _|_ q ^ £ 3^ 



6 ds 

We define the mass-like parameter \i given by 



(62) 



/! := m + 



ces 2 



ftabJ cJ 



a rich 



4 



IP T? ab 
ftnhft 



2<T 2 



B„kB 



ah 



(63) 
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This mass is indeed approximately conserved, 



Notice that p a 



ds 



mu„ 



+ O(e 6 ). 



(64) 



C(e 2 ) and © lead to to 2 = 



m + 0{e ). This implies that the Tulczyjew and the 
Frenkel condition are actually equivalent within this ap- 
proximation. 

One can guess ([63]) by realizing that - because 
of reparametrization invariance - the Lagrangian (or 
Routhian Rm) must be equal to 



R M = u a p a =m = m + £>(e 3 ). 



(65) 



In fact, (JB3J) is identical to dDTJ with S ab u b m S ab p b = 
and u 2 — 1 inserted. 



C. Dimensionless parameters and realistic values 

Before we work out the effective potential for the 
quadrupole case, it is useful to introduce dimensionless 
variables 3 , 



r 



J 



ai :-- 



J 



Mm 
a 

M' 



a 2 



E 



S 



(66) 
(67) 



Astrophysically reasonable values for the dimensionless 
spin variables are given by 



H<i. 



Notice that 



Mr 



qa 2 + 0{e 3 ), with q 



(68) 



(69) 



so spin effects are strongly suppressed for extreme mass 
ratios ?<1. 

The tidal deformation parameters \i 2 and a 2 are usu- 
ally made dimensionless with the help of the area radius 
R of the object, see [H (48) and (72)], i.e., 



k 2 := 



3GM2 
2i? 5 



.72 



48Gct 2 



(70) 



Realistic values for neutron stars are k 2 ~ 0.1 and j 2 ~ 
—0.02, while for black holes fi 2 ~ 0, and a 2 ~ 0, see 
42L l45l . For tidal deformations of black holes see also 
46H48j and references therein. An estimate for white 
dwarfs is k 2 ~ 0.01 [49] . For convenience we also define 
a dimensionless radius R, 



(71) 



3 Note that in our units c = 1, so [G] = m/kg, see also section IFl 
for an overview. 



which is just the inverse of the compactness of the test 
body. Finally, we introduce 



C 



ES 2 



(72) 



which for neutron stars is of the order Ces 2 ~ 5 [501 ] and 
for black holes it holds Ces 2 = 1- Now we are ready to 
express the effective potential in terms of dimensionless 
quantities only. 



D. Effective potential (quadrupole case) 



Let us recall how the effective potential in the pole- 
dipole case was derived. Under the assumption of equa- 
torial orbits (f5Tj) and aligned spin (f52"j). the six indepen- 
dent equations ([55]) are solved for the six variables {/>', 
p r , p^, S tr , S^, S^}. All components of p a and S ab 
can then be expressed in terms of {r, E, J, S, a, M, 
to}. Besides the rather technical issue of proofing the 
existence of equatorial orbits for our quadrupole model 
(see appendix [Cj, there are no changes in this part of the 
derivation. The spin supplementary condition in (|33[) is 
still valid/chosen, the second and third equation in (|3"3")l 
arc just the definitions of spin leng th S and dynamical 
mass m, and it was shown in [22], p. 210] that energy 
E and total angular momentum J of the test body are 
given by the same expressions as in the pole-dipole case, 
even if generic multipole corrections are included. The 
bottom line is that the solutions for p a and S ab found 
in sec. IIIIBI are still valid. The effective potential was 
defined as value of E for which p r vanishes. Notice that 
p r = implies that the orbit either has a turning point 
or is circular, see appendix [Cl 

The most important application of the effective poten- 
tial is to find circular orbit solutions. These are given by 
a minimum of the effective potential in the radial coordi- 
nate r. We therefore need to work out the r-dependence 
of the quantities entering the effective potential when our 
quadrupole model is used. The spin length S is still con- 
stant, so the only correction is coming from the dynam- 
ical mass to. Indeed, to is not a constant any more and 
given by (|63[) . More generally, we conclude that whenever 
the supplementary condition S ab pb = is used, and the 
spin length is constant, the quadrupole and higher mul- 
tipoles enter the effective potential via a r-dependence of 
the dynamical mass to only. This is similar to canoni- 
cal theories for self-gravitating bodies, in which just the 
mass-shell constraint is modified, see [5l], eq. (5.28)]. 



Inserting the solutions for p a and S ab from sec. MI Bl 
into (1631) and writing the result in terms of dimensionless 
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quantities, we obtain 



E. Tidal Disruption 



m 



C 



ES 2 Q a 2 



2f 3 
k 2 g 4 R 5 

4f 6 
-0(e 3 ). 



J - a\E 



J — a 1 E 



J - a 1 E 



J - a x E 



J — a\E 



(73) 



Consistent with our approximation we have replaced m 
by fi on the right hand side and neglected higher orders of 
S. Similarly one may replace E by the e°-order solution 
for the effective potential. Otherwise the equation p r = 
with 



aE 2 -2/3-E + 1 —~S= r 



(74) 



would contain higher orders of E, and it would be difficult 
to solve for the effective potential analytically. If one 
treats E in ([73j) pcrturbatively, then one still has only 
two solutions, 



When tidal forces become too large the test body can 
be disrupted. This limits the effects of tidal deformation 
on the effective potential. Equating the tidal force and 
the self-gravitational force of a non-rotating test body on 
its surface in Newtonian theory leads to 



™R<JL 
Gr 3 ~ i? 2 ' 



(82) 



which provides a rough estimate for the orbital separa- 
tion r at which tidal disruption may become relevant. In 
dimensionless variables this reads 



> 



2 A3 



2q z R 



(83) 



On the other hand, an estimate for ([73"]) by its leading 
post-Newtonian contribution, cf. sec. IV Dl reads (for the 
non-rotating case, &2 — 0) 



m _ k 2 q 4 R 5 



(84) 



which limits the possible difference of the masses due to 
tidal disruption as 



, 8 m 2 (3 2 — a"f J 2 fi J 
U± = ±\ r + o — + ( 75 



/' 



m ^ fc> 1 



to p r — 0. Here U± — U±/ ' [i. For completeness, we also 
give the auxiliary variables in (|T5[) in terms of dimension- 
less quantities, 



1 - 

A f a x 

'^2 



l\ qaia 2 



1 = 

S = 

a = 

A = 



A 

-2 

f 



r 



qa-2 
f 

1 + 

qa-2 
f 

2 



fO(e 3 ) 
qaia 2 



Oi 

r 



qa2 



0(e 3 ), 



5) -1 + ^ 



f 2 - 2r + a 2 . 



(76) 

(77) 
(78) 
(79) 

(80) 
(81) 



Notice that m appears here only within the S'-dependent 
terms, where it can simply be substituted by the con- 
served mass-like parameter fi within the used approxi- 
mation. 



Ii 



4 R 



(85) 



This can be understood as follows. The dynamical mass 
m also includes the tidal interaction energy, which (for 
k 2 > 0) is negative, thus the dynamical mass is reduced. 
This can be interpreted as a reduction of the gravita- 
tional potential energy (with respect to the "external" 
gravitational field) of the object due to its tidal defor- 
mation. If the tidal interaction energy reaches the order 
of fc2/(4i?), then the object is starting to get disrupted. 
This in turn limits the possible reduction of m due to 
tidal effects. Implications are discussed in the following. 

Notice that the estimate in the present section is just 
Newtonian and thus might not be accurate for very com- 
pact objects like neutron stars; see, e.g., [52[ for this case. 



V. BINDING ENERGY 

In this section we explain how the gauge invariant 
relation between binding energy and total angular mo- 
mentum can be obtained from the effective potential (for 
circular orbits and aligned spin). The various spin and 
quadrupole contributions to this relation are separated 
and compared against each other. 

The (conservative) self-force correction to the binding 
energy recently derived in (38j is also included in our plots 
as a further reference point. Results from full numerical 
simulations can be found in 1371. 
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A. Definition 

The effective potential U+ takes the value of the con- 
stant energy E of the test body in the case p r = 0, e.g., 
for circular orbits (j5"3"]) . We therefore define the binding 
energy as 

e := U+ - 1, (86) 
under the condition that the orbit is circular (|53[) . or 



dejf, J) 
df 



= 0. 



(87) 



This condition allows one to solve for the radial coordi- 
nate f, which can subsequently be eliminated to arrive 
at the binding energy e( J) as a function of the total an- 
gular momentum J. This relation e(J) is actually gauge 
invariant, which in the present context can be under- 
stood easily. That is, both the energy E and the total 
angular momentum J are scalars defined in a covariant 
manner based on the Killing vectors of the background 
Kerr spacetime, and circular orbits are a gauge indepen- 
dent concept. 

Given the square root, and the high powers of f ap- 
pearing in U + , it is necessary to solve (|8"7]l numerically. 
However, in order to better separate the different contri- 
butions to e(J) and analyze their scaling behavior in q, 
R, etc., we will expand (1551) and (|57|) in the multipole 
counting parameter e in the following. This allows one 
to solve (|87[) for r order-by-order in e analytically if the 
background is Schwarzschild, i.e. a\ = 0. We will first 
consider the Schwarzschild case in the next section and 
come back to a Kerr spacetime in sec. IV CI But in gen- 
eral it is expected to be more accurate to not expand the 
square root appearing within the effective potential, as 
this better reflects the non-linear aspect of gravitational 
interaction. This is of particular importance if one ex- 
trapolates to comparable masses q ~ 1, but this is beyond 
the scope of the present work. 



B. Corrections in Schwarzschild spacetime 

1. Spin Effects 

There is a subtlety that needs to be discussed when 
looking at effects of the test body's spin a?,. Namely, the 
total angular momentum also contains the spin of the test 
body. In canonical formulations, the total angular mo- 
mentum is just the sum of the canonical orbital angular 
momentum l c and the spin of the object, i.e., 

lc ■■= j^(J - S) = J - qa 2 . (88) 

The reason is that the total angular momentum generates 
rotations of the whole system on the phase space, which 
directly decompose into a rotation of the orbital variables 



(generated by l c ) and a rotation of the object (generated 
by the spin), see, e.g., [HI, sec. 4.1.2]. If one now considers 
the Newtonian energy of point-masses 



(89) 



4_i, 

2f 2 f ' 

and rewrites it in terms of the total angular momentum 



j2 i 



a 2 Jq 



OU 2 



(90) 



an apparent spin-dependence of the binding energy 
arises. However, if one compares spinning and non- 
spinning systems at the same total angular momentum 
J, one is comparing them for different orbital angular 
momentum, i.e., one is comparing them in different or- 
bital configurations. The difference in their energies is 
just due to this change of the orbit, and not due to an 
actual spin interaction. There is nothing wrong here, this 
just shows that comparing e(J) can be misleading. We 
will therefore use e(l c ) for comparisons here, as for the 
same value of l c the system is approximately in the same 
orbital configuration. 

In order to better separate the different contributions 
to the binding energy, we make a series expansions of the 
form 



r = r + e r\ 



e 2 f 2 



(91) 



This, together with J = l c + qa 2 , is inserted into ([87]) 
and the whole expression is expanded in e. Then one can 
solve for f order-by-order in e, i.e., the e°-part of (1ST]) is 
solved for r , the e 1 -part of ([87]) is solved for ?i, and so 
on. The solution for f is then plugged into e(f, l c ) leading 
to e(l c ), which is expanded again, 



e(/ c ) = e (/ c ) + e e x {l c ) + e 2 e 2 (l c ) + 



(92) 



The correction e\(l c ) does not contain quadrupole contri- 
butions and is linear in the spin a 2 and mass ratio q. In 
the figs. 2(a) and 2(b) we have plotted the "normalized" 
relative difference to the leading order eo(l c ), 



logi 



ei(Zc) 



e (l c ) 



logio 1 902 1 



(93) 



which does not depend on a 2 and q (this corresponds to 
setting a,2 = 1 = q). The correction e 2 (l c ) can be split as 

c 



e 2 {l c ) 



('=)• 



(U + 4 2 (0 + 4 2 (Z c ), (94) 



into tidal quadrupole contributions e 2 {l c ) and ej 2 (Z c ) 
proportional to k 2 and j 2 , respectively, a quadratic-in- 
spin quadrupole e 2 Bs2 (l c ) proportional to Ces 2 > an d a 
remaining term ef (l c ) which is quadratic in spin. 
The spin parts are again normalized, 



log 



10 



r '2 



(lc 



eo(lc 



log 10 \q 2 al\ 



log 



10 



'(lc) 



eo(Z c ) 



lot 



10 



\n 2-21 
\^ES 2< 1 a 2 \ 



(95) 
(96) 
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and plotted in figs. 2(a) and 2(b) - tidal effects will be dis- 



cussed in the next section. These plots also show recent 
results for the conservative part of the self-force normal- 
ized as 



lot 



10 



eTBB Qc) - eo(l c 



e {l c ) 



log 



10 



(97) 



where eTBB A) ^ s ^ ne binding energy given in a paramet- 
ric form in [38|], see also app.lElfor an explicit expression. 
The difference exBB('c) — eo(l c ) changes its sign at a cer- 
tain value of l c , which appears in the plots as a pole. It 
should be stressed that our results do not include contri- 
butions that arise from a perturbation of the background 
spacetime due to the spin or quadrupole of the test body. 

Notice that the transition from the normalized curves 
to realistic curves in figs. 2(a) and |2(b)| is straigtforward. 
If one considers a mass ratio of q = 10 -2 , then the self- 
force curve 0971) is just shifted down by log 10 1 = — 2. 
Similarly, for a rotating star - with &2 = 10~ 05 s» 0.32 
and Ces 2 — 10 ' 5 3.2 - the correction linear in spin 
(|9"5]l is shifted down by —2.5, the curve 09"5j) is shifted 
down by —5, and the Ces 2 -quadrupole curve (l96|) is 
shifted down by —4.5 (still for q = 10 -2 ). For a black 
hole one may even assume hi ~ 1, furthermore we have 
Ces 2 = 1- The latter case is illustrated in fig. 2(c) No- 



tice that the dominant scaling is coming from the mass 
ratio q. The linear-in-spin correction and the self-force 
scale with the same power of q. One can also easily ob- 
tain the sign of the correction Ae to the binding energy, 
it holds 



sgn (ei) = sgn (a 2 ) , sgn (ef J 
sgn (e^ Es2 ^ = -sgn(C E s 2 ) , 



-1, 



(98) 
(99) 



while the self-force correction exBB — eo has a positive 
sign for small l c and a negative sign for large l c . 

The plots confirm that the multipole approximation 
introduced in the present paper is well justified, as the 
quadrupole effects are always at least an order of mag- 
nitude weaker than the linear spin effects - within the 
context of the quadrupole model adopted in this work. 
Also notice that the quadrupole effects scale with a higher 
power of the mass ratio q than the linear spin corrections. 

To further illustrate the difference between the param- 
eters l c and J, we also include a plot of the binding energy 
e( J) in fig. m 



2. Tidal Effects 

For the tidal contributions we also plot normalized 
quantities, 



logio 
iogio 



eo(l c ) 



e (l c 



logio 
!ogi 



k 2 q 4 R 5 



(100) 
(101) 
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(a) Normalized curves for spin and self-force corrections. 
-0.5 




4.5 5.0 

lc 

(b) The same as|(a)| but zoomed in at small l c 



6.0 




5.5 



6.0 



(c) The same as above, but the normalized curves are shifted 
down to a mass ratio q = 10~ 2 for 8,2 = 1 and C ES 2 = 1, which 
can occur for a black hole. 

FIG. 2: Plots illustrating spin and self-force corrections to the 
binding energy of a m ulti polar object in Schwarzschild space- 
and (b) the corrections are normalized 



time ai = 0. In ( 



by subtracting their dependence on q, 0,2, and C ES 2, see (I97|) . 
flgg), (USD, and 055} [belonging to curves (*), (**), (***), and 
(****), respectively]. 
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FIG. 3: Binding energy in terms of j. This plot makes clear 
that one has to be careful when in comes to the parameter- 
ization of the strength of different corrections. In contrast 
to the Z c -parameterization, effects from the internal structure 
- in this case the spin - of the test body, appear more pro- 
nounced in the J-parametrization. The quadrupole curve is 
not affected. 



see figs. 4(a) and |4(b)[ Notice that the change in the 
effective potential due to tidal deformations, and thus 
their backreaction to the orbit, is proportional to the dif- 
ference between dynamical mass m and constant mass 
parameter fi. The scaling in q and R then immediately 
follows from (l73l) . This suggests that tidal effects are 
most relevant for big (non-compact) objects due to the 
proportionality to R . However, due to tidal disruption 
the maximal possible tidal backreaction is proportional 
to l/R, cf. eq. ([85]). Therefore, in order to produce a 
strong tidal backreaction, it is in fact important to con- 
sider objects which are very compact (like white dwarfs 
or neutron stars) , such that their own gravitational held 
can still hold the object together when the tidal forces 
become strong. 

For example, for a white dwarf one can take ki = 0.01 
and R = 10 4 . For a mass ratio of q = 10~ 4 the k^- 
quadrupole curve is then shifted up by +2, while the 
self-force curve is shifted down by —4. Tidal disruption 
can be expected when the fc2-curve reaches 



-6.6, 



(102) 



see (|85j) . Thus, the tidal backreaction can reach the same 
order as the self-force before the object is disrupted. For 
neutron stars one has R ~ 5 and ki = 0.1, so disruption 
only happens at about —2.3. However, this Newtonian 
estimate might not be very accurate for neutron stars). 
On the other hand, the /c2-curve for a neutron star is 
shifted further down by —1.5 even for a mass ratio q = 
10 _1 (where the used approximation is not accurate any 
more). For neutron stars the tidal backreaction is most 
relevant for comparable masses, as it can become very 



self-force (*) 
&2-quadrupole (**) 
y'2-quadrupole (***) 





-2 




-4 




-6 






< 


-8 


o 




an 




_o 


-10 




-12 




-14 



10 



15 



20 



I, 



(a) Normalized curves for tidal and self-force corrections. 
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c 




3.5 4.0 4.5 5.0 5.5 

l c 

(b) The same as | (a) | but zoomed in at small l c 



6.0 




lc 

(c) Similar as above, but the normalized curves are shifted down 
to a mass ratio q = 10 — 4 for &2 = 0.01 and R = 10 4 , which can 
occur for a white dwarf. Tidal disruption can be expected when 
the fc2-curve reaches —6.6. 

FIG. 4: Plots illustrating tidal and self-force corrections to 
the binding energy of a multipolar object in Schwarzschild 
spacetime fii = 0. 



In (a) and (b) 



the corrections are nor- 
malized by subtracting their dependence on q, R, hi , and ji , 
see (§7}, (fTUuT) . and (THJTj) [belonging to curves (*), (**), and 
(***), respectively]. 
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(a) Linear in spin ~ I12, eq. d 93 1 ) 
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(b) Quadratic in spin ~ S|, cq. {95} 



FIG. 5: Pole-Dipole corrections in Kerr spacetime. 



strong in this case. The signs of the corrections are given 
by 



sgn ( e 2 2 ) = - s g n , sgn 



C. Corrections in Kerr spacetime 

For Kerr spacetime, ai ^ 0, one can still separate the 
various contributions to the binding energy. This is com- 
pletely analogous to the case oi = 0, except that we are 
not solving for tq analytically. Still one can write the 
zeroth-order approximation of ([57]) as [53, eq. (2.13)] 



l,= 



r| - 201^ + 0^ 
/ru"^ ?l + 2a 1 y/¥r] - 3f 



(104) 



which is most conveniently inverted numerically. 

Plots for the linear and quadratic spin contributions for 
a pole-dipole particle are shown in fig. [5J and quadrupole 
contributions can be found in fig. [5] For negative a\ the 
contributions are stronger than the Schwarzschild ones 
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(a) C ES 2 -quadrupole proportional to Sj, eq. H96I I 
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(b) Tidal A;2-quadrupole, eq. (1 100 1) 
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(c) Tidal j'2-quadrupole, eq. 110H 
FIG. 6: Quadrupole corrections in Kerr spacetime. 



for the same value of l c - However, the last stable circu- 
lar orbit is also reached at greater values of l c , so finally 
larger corrections are possible in the Schwarzschild case 
(at values of l c where the Kerr orbits are already unsta- 
ble). For positive d\ the opposite is true: the curve is 
below the Schwarzschild one, but the last stable circular 
orbit is located at lower l c . For the case a± = 1 the last 
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stable orbit reaches the horizon of the black hole. There- 
fore the test body can orbit much closer to the horizon, 
where the held is very strong. This finally allows the same 
contributions to become one or two - or even more for 
tidal effects - orders of magnitude stronger than in the 
Schwarzschild case if Si > 0. This can make quadrupole 
effects for neutron stars interesting again even beyond 
the comparable mass case. 

D. Comparison with post-Newtonian Hamiltonians 

We will now compare our result for the binding en- 
ergy e(J) with post-Newtonian Hamiltonians. We first 
calculate the post-Newtonian expansion of our result for 
e(J) (which is valid to all post-Newtonian orders). Then 
again we calculate e( J) from post-Newtonian Hamilto- 
nians (valid for generic mass ratio q) and expand in 
the mass ratio q. Both should lead to the same re- 
sult, as the regimes of approximation overlap and the 
relation e(J) is gauge invariant. We include in our 
comparison all (conservative) post-Newtonian Hamilto- 
nians in the Arnowitt-Deser-Misner transverse-traceless 
(ADMTT) gauge known to date. 

We start with a definition of the post-Newtonian ex- 



pansion. The post-Newtonian approximation is a weak 
field (f 3> 1) and slow motion (l c ^> 1 or J 3> 1) approx- 
imation. For bound orbits the Newtonian virial theorem 
establishes a relation J 2 ~ f. It is convenient to intro- 
duce a post-Newtonian book-keeping parameter epN and 
associate power counting rules to all variables, 

l c = 0{e- p x N ), J = 0{ep 1 N ) 1 r = 0(ep%). (105) 

Notice that Si, a 2 , C E s 2 , k 2 , J2, and R are in fact fur- 
ther expansion variables that are not related to the post- 
Newtonian approximation, which suggests to count them 
as 0{e PN ). This approach will indeed be followed here 
with the only exception that we formally assume 

R = 0{e- p x N ), (106) 

This choice is of course quite arbitrary and is only made 
because otherwise the tidal interactions would formally 
appear at very high post-Newtonian order only (which 
somewhat underestimates their effect and would require 
us to expand to very high order for the comparison). 

Using these power counting rules it is straightforward 
to expand e(f, J) in epjv, 
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32f 3 128f 4 256f 5 2r 
3a 2 a 2 J 2 g 2 3aia 2 g 9ciia 2 J 2 q 15aia 2 J 4 q 



a 2 J 2 q 2 



+ - 



8r 

- OJ02 
J(? 4 i? 5 

4f 3 



13) 



2f 
16f 2 



a 2 jV 



(9C BS2 + 17) 



2r + 


8f 2 


7a 2 J 6 g 2 


axJq^R 5 


32f 3 


h 2f 3 




- -1^ 


2Jf 2 


4Jr 2 



4f 

"2 2 

a2 9 



(C 



4f 2 16f 3 
3) 



ES 



(12fc 2 +j 2 ) 



(12fc 2 +. 72 ) 



(10fc 2 +j 2 )- 



J 3 q 4 R 5 
8r 4 



■(15fc 2 +j 2 ) 



(107) 



Next we solve the condition defining circular orbits (|87[) for f (order- by-order in epjv) and insert the result into (|107l) . 
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This leads to the post-Newtonian expanded gauge-invariant relation 

1 



/ft = \_ «2g 1 

2J 2 J 3 J 4 



9 3*2 2 

8 ~ r* q 



1 

J 7 
1 

1 

J9 
1 

J 1 " 



1 

2„2 



I8&1 — 185.2(7 — 3a 1 a2<7 + 2Aa\a 2 q — k 2 q 



2ai — 30,2^ 



81 1 1 

- — - -a\ + 9aia 2 g - -a|g 2 (2C BS 2 + 15) 



J 5 



3861 89 



128 



183 21 n n n 21 n n . l 



4J 5 



162ai + lQaf - 135a 2 o - 151a 2 a 2 g + \2a x alq\C BS 2 + 20) 



•(78fc 2 +j 2 ) 



53703 5643 . 3„ 4 7623 

a, a, H 

256 16 1 2 1 8 



aia 2 g + 84a^a 2 g 



~a 2 2 q 2 (182C ES2 + 907) - ^a 2 ^^ + 125 ) + ^-^!( 8 4fc 2 + j 2 ) 



1 

J 11 



1512ai + 386aJ - 1134a 2 g - 2619a^a 2 g - 15a?a 2 g + 6aia|g 2 (37C ElS 2 + 453) + 390a?a 2 q 2 



47 



1 



g 4 i? 5 / 224 7, 



(108) 



Now we establish the connection to the (conserva- 
tive) post- Newtonian Hamiltonian H within the ADMTT 
gauge. We include the Hamiltonians 



H = H* 



Hi 



PN 



2PN 



H. 



3PN 



H 



so 



i ttSO 
+-"NNLO 



H 



s 1 s 2 

LO 



H 



SiS 2 
NLO 



LO 

S1S2 
NNLO 
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SO 
NLO 



H 



LO 



+-^NLO + ^LO 



H 



LO 



^LO 



(109) 



which we will list in the following in the center-of-mass 
frame for aligned spins together with the corresponding 
literature. Further, we already expand the Hamiltonians 
in the mass ratio q to the order needed for the compari- 
son. 

The Newtonian (N), first post-Newtonian (1PN), and 
second post- Newtonian (2PN) Hamiltonians are given by 



-HlPN 
-ff2PN 



1 

1 

2ft 



1 - 



1 - 



C L 

1 



1 - 



/2 
l c 

2f c 

f c 4f 2 
10Z|_ 5^ 
f r 2ft 



I 6 

Aft 



(110) 
(111) 
(112) 



see, e.g., [5J-|56j, and references therein. Here r c denotes 
the ADMTT-gauge (canonical) radial coordinate. The 
third post-Newtonian level was first tackled in [53] • But 
the result contained two "ambiguity" parameters , which 
were subsequently determined [58|, |59| , see also [60] • (Di- 
mensional regularization must be used to avoid such am- 
biguities (6l|.) An alternative derivation of the equations 
of motion at the third post-Newtonian order can also be 
found in [62|, |63| . The full result in the test body limit 



reads 

#3PN 



1 

8ft 



mi 

fr 



27lt 
2ft 



2ft 



5/ 



8 1 



16ft 



(113) 



The leading order (LO) spin-orbit (SO) and SiS 2 Hamil- 
tonians are given by 



rrSO 



TTS1S2 

n LO 



2a\l c 3qa 2 l c 



2ft 



qa\a 2 



(114) 
(115) 



see [64-67] . These references also contain the S\S\ inter- 
action potential for the black hole case Ces 2 = 1- The 
parameter Ces 2 was introduced in [68j . The sum of Si Si 
and S2S2 interactions can be written as 



H s - 



2f 



Ces 2< 1 ^2 
2ft 



(116) 



The next-to-leading order (NLO) spin-orbit Hamiltonian 
was derived within the ADMTT gauge in [69l.l70l]. 



rrSO _ 
-"NLO — 



6ail c hqa 2 l c 



r 2 

8ft 



(117) 



Other results on this interaction can be found in [ 7l |-[76l . 
The next-to- leading order Si S2 Hamiltonian [70, |77| is 
given by 



H 



S1S2 _ 6gaia 2 



NLO 



I 2 

1- -2- 

4r fi 



(118) 
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see also [7J, |7J, |79( . The next-to-leading order Si Si and 
S 2 S 2 interaction Hamiltonians reads [80l. l8ll 



-"NLO 



2 ~ 2 
f«2 



2rf 



(2C ES 2 + 1) - 
5 



HC 



ES 2 



4f R 



3£ 
2f, 



(119) 



see also [23, [26|, [76|, [83. Recently even next-to- next- 
to-leading order (NNLO) spin interaction Hamiltonians 
were calculated, namely the spin-orbit one in the test- 
spin limit [36j (see [83| for the complete result) 



TT 



SO 

NNLO 



21di^ c 3ga 2 Z c 



2f 5 c 8f5 
and also the Si S2 Hamiltonian [8J, [8 



(120) 



TTS1S2 

"NNLO 



4f5 



4; 2 

rv 



2f? 



(121) 



The Hamiltonians cubic in the spins derived in [86|, |87J 
are only valid for binary black holes (i.e., Ces 2 = 1) an d 
can be summarized as 



#lo = ^5 (Si + <?a 2 ) 2 (4ai + qa 2 ), 



(122) 



see also [3a]. 4 The leading order tidal interaction Hamil- 
tonians read 1281 



LO 



k 2 q 4 R 5 



Tin 



hq A R 5 l 2 
4f? 



(123) 



Hig her order tidal interactions are considered in [28|, [29|, 
l88j , but results therein can not immediately be included 
in the present comparison as they are not given in the 
form of Hamiltonians within the ADMTT gauge. 

Finally, we solve the condition defining circular orbits 



dH{f c ,l c 

8f n 



0, 



(124) 



order-by-order in epN for the ADMTT radial coordinate 
f c and eliminate f c from the Hamiltonian. We find that 
the result H{l c ) agrees with e(J) given by f)108|) for all 
the Hamiltonians shown in the present section (taking 
into account that l c = J — ga 2 ). 



4 Also Hamiltonians of quartic order in spin are given in |8a . |87|1 . 
However, complete agreement with the results of the present pa- 
per could not be found yet. The deviation indicates that the S 4 
Hamiltonian discussed in [8^1 does not vanish, in contrast to the 
conclusion therein. Also notice that a misprint at quartic order 
in spin was corrected in the arXiv version of [gij recently. 



VI. CONCLUSIONS 

In this work we have investigated the influence of the 
internal structure of test bodies on their motion within 
the context of a multipolar approximation scheme. Cor- 
rections arising from the spin (dipole) and quadrupole 
moment were worked out explicitly for aequatorial orbits, 
with aligned spin, in Kerr spacetime. In particular our 
explicit model for the quadrupole, which also allows for 
tidal deformations of the test body, goes beyond previous 
investigations in the literature in a multipolar context. 

Our comparison with recent numerical results for 
structureless bodies - which take into account the self- 
force - makes clear, that the corrections arising from the 
extendedness of the body can play a role in the descrip- 
tion of the motion and should be carefully dealt with. A 
final statement regarding the magnitude of effects coming 
from different corrections (internal structure or self-force) 
depends on many factors. The figures in the present work 
immediately identify the relevant contributions needed to 
achieve a specific accuracy in the binding energy. As we 
have pointed out in the context of the binding energy, one 
should be careful when it comes to the parametrization 
of possible contributions in terms of different variables, 
cf. the corresponding discussion regarding parameters l c 
and J. One should also stress at this point, that only 
the conservative parts of the self-force in Schwarzs child 
spacetime - which were the only ones readily available in 
the literature - were taken into account in our compari- 
son. 

Our work clearly indicates the necessity of future in- 
depth comparisons of different approximation schemes. 
In particular, as soon as results for the self-force in Kerr 
spacetime become available they should be compared to 
our findings here. Another interesting open question re- 
gards the possibility of extrapolating our results to less 
extreme mass ratios, i.e., intermediate mass ratios or even 
comparable masses, in the future. Whether such an ex- 
tension of one of the existing multipolar approximation 
schemes can be consistently worked out is open for de- 
bate. A promising approach is to substitute masses or 
other parameters in a certain way, which at least for the 
self-force works astonishingly well [8!| and can also be 
understood as a change of expansion parameters [38| . A 
comparison with post-Newtonian results (not expanded 
in the mass ratio) can serve as a guide to identify proper 
variable replacements. The post-Newtonian expansion of 
the gauge invariant relation e(J) can further be used to 
match coefficients of an effective action with higher order 
spin couplings in the future. For example, if we would 
have kept the constant Ces 2 f° r both objects in the post- 
Newtonian Hamiltonian Hf or H§ LO (which can be de- 
rived directly from the effective action [25|, [5l[ ) , then a 
comparison with the result from the present paper would 
show that one has to set Ces 2 — 1 f° r the central black 
hole. The relation e( J) is also very useful to check results 
of post-Newtonian or post-Minkowskian approximations. 

The abovementioned extensions of the approximation 
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method are required for the creation of gravitational wave 
template banks in the comparable mass case. Although 
numerical simulations are optimally suited for the late 
inspiral phase of such a binary, they currently are not 
able to cover the whole parameter space in an acceptable 
timeframe if both objects are spinning. A synergy of 
numeric and analytic methods is needed. 
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Appendix A: Misprints within the effective potential 

There seem to be minor errors or misprints in the lit- 
erature on the effective potential for the pole-dipole case 
[3ll - l34l |. We will summarize here what misprints needs 
correction in order to achieve agreement with the result 
in the present paper, equations (Pd| - (I5(I1) . and (1521 . 

Equation [3l], (7a)] must read 



= A 



2sa 



3 + — 



(Al) 



where the dots are an abbreviation for correct terms in 
this section. In [34j there are two sign errors in the Ap- 
pendix, namely 

S 2 AM 2 , 
J 



7 

k; 



Bi 



Mr 



(A2) 
(A3) 



In [33j, eq. (2.24)] the expression for Z must read 

2 



Z = - 



(A4) 



Finally, we find full agreement with [321 ] except for an 
overall sign of the spin S. We will argue in the next 
section how to correctly identify the co-rotating and 
counter-rotating cases. 



Appendix B: On the orientation of the spin 



tation depends on the sign choice in (l34l) . which in turn 
depends on conventions for the spin tensor and the sig- 
nature of spacetime. 

The simplest way to identify the spin orientation for 
the sake of the present paper is via the angular momen- 
tum J defined by 



J 



+ 



pt,rt 



rPt) 



S 

2mr 



(Bl) 



In the weak field (large r) and slow motion limit we have 



9^t,r ~ 0, g^ >r s 
Therefore it holds 



-2r, 



-mr <p, and pt ~ m . 



J&L + S, 



(B2) 



where the orbital angular momentum L = mr 2 <t> is 
aligned to d z if L > 0. This shows that for S > the 
spin is aligned to d z . 

More generally, one may define the angular momentum 
J in the weak field and slow motion limit via the energy- 
momentum tensor as 



r 



tiki 



d 6 xx k T 



krplO 



(B3) 



which is still applicable if the spacetime possesses no ro- 
tational symmetry. From (jSJ we obtain 



1 



T' u « mY l S(x l - Y l ) - -OkiS^Six 1 ~ Y 1 )), (B4) 



and J3S leads to S l w he ik iS kl . Therefore it holds 



J' 



S' 



(B5) 



where L % := meikiY k Y l is the usual Newtonian orbital 
angular momentum vector. Notice that for usual spheri- 
cal coordinates dg points in the opposite direction as d z 
in the equatorial plane. Therefore S < corresponds to 
a spin aligned to d z . 



Appendix C: Existence of Equatorial and Circular 
Orbits 



The conditions for equatorial orbits (|3T|) and aligned 
spin (I3"2l must be preserved under the time evolution, 
i.e., 



= u a = 0, f = 0, S aS = 



(CI) 



Similarly, for the interpretation of the effective potential 
and the existence of circular orbits, it is important that 
under the condition p r = it follows that 



0. 



(C2) 



Whether the spatial components of the spin vector S a 
allow a straightforward determination of the spin orien- 



or in words, if p r = then there is either a turning point 
of the orbit or the orbit is circular (in both cases u r = 0). 
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In order to prove these statements one can use the re- 
lation (I16[) and the equations of motion. It is further 
beneficial to use the coordinate time as the worldline pa- 
rameter instead of the proper time, i.e., u l = 1 (which is 
why we formulated our model (|54[) in a reparametrization 
invariant manner) . In the pole-dipole case the statements 
above where shown in, e.g., [U, [34j]. This calculation 
must be repeated for the quadrupole model (|54l) now. 
As the details of such a calculation do not provide any 
physical insight, we simply relied on a brute force cal- 



culation using Mathematica 90] together with the free 
xTensor and xCoba packages 



91 



We found that the 



statements generalize to the quadrupole case. 



Appendix D: Quadrupole action 

In this section we show the extension of the point-mass 
Lagrangian that inspired the quadrupole model used in 
the present paper. It is actually more convenient to work 
with the Legendre transformation of the Lagrangian in 
the angular velocity, so we are technically considering a 
Routhian Rm here. Combining models for quadrupole 
deformations induced by spin from [2^, eqs. (1) and (16)] 
and for adiabatic tidal quadrupole deformations given by 
[H eq. (19)] or by [H eq. (5)] leads to 

Rm = ^ ^Ba^UcS^ - ^$=E ab S a c S cb 



/'-2 



4(Vu 2 ) 3 



IP J? ab 



2cr 2 „ ,-,„,, 



3(vV) 



B ab B ab , (Dl) 



where u 2 := u a u a . Notice that some signs changed due 
to the adoption to our conventions. The parameters /i, 
c-es 2 i M2j and 02 are assumed to be constant. One may 
write fi = mo + jjS 2 + . . . , where / can be interpreted as 
a moment of inertia and mo as an irreducible mass, see, 
e.g., 0, eq. (3.28)]. 

The equations of motion for Lagrangians of the type 
used above have been worked out in a general fashion 



already in [4(| and were found to be of the form used in 
the present paper (see also sec. 5.2]). In particular, 
[iol eq. (19)] provides a formula for the quadrupole J abcd : 



jabed = g 



dR 



M 



dR a 



bed 



(D2) 



It is straightforward to derive J abcd from this formula. 
But the equations of motion belonging to the Rm, 
shown in (|D1[) . preserve the spin supplementary condi- 
tion S ab pb — only approximately (in the sense of the 
multipole approximation introduced in sec. HVBj) . while 
here we are enforcing this condition exactly. Therefore we 
made some minor changes in the result for J abcd } which 
are, however, in accordance with the used approxima- 
tion scheme. That is, we replaced u a by p a /m and in- 
troduced the overall factor m/m in (f54"|) for the sake of 
reparametrization invariance (which was ensured by fac- 
tors of y/u? in the original expression) . 

Appendix E: Explicit expressions for the binding 
energy 

In this appendix we provide explicit expressions for 
the expanded binding energy if a\ = 0. Although these 
can be easily derived from the more compact effective 
potential given in the main text using computer algebra, 
we display them here for the sake of completeness. The 
zeroth order solution to (l87l) reads 



(El) 




and the corresponding binding energy is given by 




(E2) 



It is straightforward to solve for higher orders in terms 
of ro and eo, resulting in 
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e(/ c ) 



eo + 



~-3/;2 , -2^1 ro;2/~ ,< \ 0-21-! J,),,*. .j,2,,. I 



f -2 

«2f2(2 + e f )(6 + f2-6fo) 



C(f - 3^ (f - 2) + [6 + (2 + 2e - 3)f ] 



^2 2 

a29 



+ qM*o - 2)^f - a (^ + fg)" a (2^ - 3/ 2 (f - 4))"^ 8(4 + 4e - 3)(-2 + f ) 2 fj° 



3Zj°(f - 2) 



12(r - 2) 2 [16 + (f - 9)f ] + (1 + e )[f (449 + f (4(31 - 3f )f - 399)) - 110] 



2(f - 3) 2 (f - 2) [24 + (f - 8)f ] - (1 + e )[f (3627 



;4 f 6 
'c'O 



6(f - 4)(f - 2) 2 [116 + f (4fo - 49)] 



+2f (f (314 + 3(f - 16)f ) - 1056)) - 2484] 

-(1 + e )[f (6200 + f (-5601 + 2f (1087 + 4f (3f - 47)))) - 2220] 
-lift 6(f ~ 2) 2 [f (300 + f (-37 + (f - 6)f )) - 504] - (1 + e )[10116 + f (f (11801 
+f (f (1027 + 6(f Q - 21)f ) - 4520)) - 17028)] 

6(f„ - 2) 2 [108 + f o (f o (105 + 2(f - 12)f ) - 192)] - (1 + e )[26892 + f (f (41277 



'c'O 



+f (r (4115 + f (27f - 520)) - 17358)) - 52056)] 

1 + eo 
ro -2 



fl{ll + fl)- 1 [2fl 3/ c 2 (f - 4)] 1 |2(f - l)f 2 + ll [(16 - 3f )fo " 18] } ( 

+ 0(e 3 ). 



m 



-(l + eoyUlt + ftr'm-^ + ft] [3/ 2 (4-fo) + 2f 2 ]- 1 ^ ■ 



(E3) 



where the mj \i contributions are given by (|73|) with J ss 
! c , r w fo, and « 1 + eo inserted. 

The self-force correction can be derived from the for- 
mulas in [2U and reads 



6TBB (^c) - e (l c ) 



1 



1 



i/l - 3a; 2 



where 



6*o = l-Jl-f. 



zsf(^o), (E4) 



(E5) 



The function zsf(x) is given in [38| by a fit to a rational 
function with 5 parameters. 



r 



different quantities appearing throughout the work are 
displayed in table U and table Ql] Note that we set c = 1 , 
the dimension of the gravitational constant then becomes 
[G] = m/kg. Table IILTI contains a list with the most im- 
portant symbols used throughout the text. Latin indices 
denote 4-dimensional indices and run from a = 0, . . . , 3, 
the signature is (+,-,-,-). The Riemann tensor R a bd c is 
defined by 



V[ a V fc ja c = -Rab/a , 



(Fl) 



Appendix F: Conventions & Symbols 

In order to fix our notation, we provide some tables 
with definitions in this appendix. The dimensions of the 



where a c is a generic vector. The volume form is given 
by riabcd = y/^ge a bcd, where e abc d is the completely anti- 
symmetric Levi-Civita symbol with 60123 = +1- Notice 



that rf 



abed 



e abcd /y/=g and e° 



1. 
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TABLE I: Dimensions of the quantities (general). 
Dimension (SI) Symbol 



Geometrical quantities 


1 


9ab, \J-g, o b , rj , e , £ 


m 


s, Y a , dx a 


m" 1 


d a , r a (, c 


m- 2 


Rabc , E a b, B a b 


Matter quantities 


1 


u a , U a , p a , H 


kg 


m, m, p a , p, R M 


kgm 


S ab , S a , S 


kgm 2 


jabcd Qa&c 


kg/m 3 


rpa.fi 


Auxiliary quantities 


1 


P a , r 


m- 4 


8(4) 


Operators 


1 


pi, W b a , p a b , X a b 


m" 1 


vv D a • „ 


TABLE II: Dimensions 


of the quantities (Kerr). 


Dimension (SI) 


Symbol 




Geometrical quantities 


1 


9, 4> 


m 


t, r, M, a 


Matter quantities 


kg 


E, Pt , p r , S**, S r f, U± 


kgm 


J, p<p, pe, S tr , Se 


Auxiliary quantities 


1 


a, P, 7, 5, a, e, C ES 2, R, U 




f, J, E, ai, a 2 , q, 32, k 2 , l c 


m 


P 


m 2 


A 


kg" 2 m~ 5 


A 


kg 2 m<"> 


,(n=3) .(4) .(5) .(7) 
J± 1 , I\ 2 , -^3 i -^4 


kg" 1 


C ES 2 


kgm 4 


P2, U2 
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TABLE III: Directory of symbols. 



Symbol Explanation 



Geometrical quantities 



9ab 

c 

x a , s 
Y a 

Rabc 1 
E a b, B ab 

M, a 



Matter quantities 
~ a 

m, m 

/•< 
P a 

s ab , s a , s 

rpab 

Rm, H 



Auxiliary quantities 

[11, &2, C ES 2, C ES 2, J2, fc 2 

R 
q 

u±, u± 



Operators 

Pb, xs, w b a 

abed abed 
if , fc 

ft, v, 

_D _«■„ 

ds 



Metric 

Determinant of the metric 

Kronecker symbol 

Killing vector 

Coordinates, proper time 

Worldine 

Connection 

Curvature 

Curvature (electric, magnetic) 
Kerr (mass, parameter) 



Velocity 

Mass (Frenkel, Tulczyjew) 
Mass-like parameter 
Generalized momentum 
Spin (tensor, vector, length) 
Energy-momentum tensor 
Routhian, Hamiltonian 



Coupling constants 
Dimensionless radius 
Mass ratio 
Effective potential 



Spatial projectors 
Permutation symbols 
(Partial, covariant) derivative 
Total derivative 
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